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Abstract 

We point out that if heavy metastable particles composing the dark matter of our 
galaxy are responsible for the ultra-high energy cosmic rays (UHECR) then the leading 
tidal stream of the Sagittarius dwarf galaxy could be detected through UHECR. The 
signal would be an anisotropy in the UHECR flux smaller than the telltale anisotropy 
towards the galactic center that would first establish unstable dark matter as the origin 
of the UHECR. 

Cosmic ray particles with energies above the Greisen-Zatsepin-Kuzmin cutoff of about 
5 X 10^^ eV have been detected by a number of independent experiments over the last decade 
The existence of these ultra-high energy cosmic rays (UHECR) presents us with a 
problem. Nucleons and photons with those energies have short attenuation lengths and could 
only come from distances of 100 Mpc or less, while plausible astrophysical sources for those 
energetic particles are much farther away. Recently, results from the HiRes Collaboration 
|3] brought the violation of the GZK cutoff into question and this issue will only be resolved 
conclusively by the Pierre Auger Observatory. 

Among the solutions proposed for the origin of the UHECR is the decay of supermassive 
relic particles which may constitute all or a fraction of the dark matter 4j. Such particles 
must be metastable, with lifetimes exceeding the age of the Universe. In this scenario the 
flux (j) of UHECR coming from a particular volume element in the Universe is proportional 
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to the density of dark matter in this volume, and inversely proportional to the square of the 
distance r between us and this particular volume element. Thus the flux is the integral of 
the column density over the solid angle 
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Consequently, in decaying dark matter models the UHECR are produced mostly in our 
galaxy and the key test is the expected anisotropy in the arrival directions caused by the 
offset of the Sun from the center of the galaxy [S]. In these models cosmic rays from 
dark matter dominate at energies above about 6x10^^ eV, and are mostly photons [6j. 
Even for protons the effects of magnetic fields on UHECR trajectories can be neglected 
at these energies jZj, thus the trajectories of these UHECR are straight lines and point 
towards the place of origin. Spherical halo models yield an anisotropy symmetric around 
the direction of the Galactic Center (GC) (while triaxial models, which will not be considered 
here for simplicity, give considerable angular deviations). The amplitude of this anisotropy 
is controlled by the halo core radius i?c, a parameter which appears in the halo density 
models, the simplest of which is the isothermal halo model 



where r' is the radius from the GC. 

The halo core radius Rc is not well determined from astrophysical means and could in fact 
be measured through the anisotropy itself. In isothermal halo models, the GC asymmetry 
amplitude, AcciJ^)) defined as the ratio of the flux in the direction of the GC and the flux 
from the galactic anti-center (AC) within an integration cone of aperture 9, depends strongly 
on Rc. For example in [S], for i?c = 5 kpc, the amplitudes are found to be (see Fig. 5 of [S]) 
Agc(5°) ~ 5 and ^4^^(80°) - 3, while for R^ = 10 kpc, Agc{5°) ^ 3 and Agc(80°) ~ 2. 

The observation of a GC-AC asymmetry would prove the decaying particle scenario for 
the origin of the UHECR. In this case, the leading tail of the Sagittarius Dwarf Galaxy 
would give rise to a different asymmetry in the UHECR fluxes. This is the new asymmetry 
we are considering here. 

The Two Micron All Sky Survey (2MASS) and the Sloan Digital Sky Survey (SDSS) 
CD] have traced the tidal stream of the Sagittarius (Sgr) dwarf galaxy more than 360° 
around the sky. This galaxy is a dwarf spheroidal of roughly IO^Mq. It is a satellite of the 
Milky Way Galaxy, located inside the Milky Way, ~12 kpc behind the Galactic Center (GC) 
and ~12 kpc below the Galactic Plane (GP) ^T]. There are two streams of matter that 
extend outwards from the main body of the Sgr galaxy and wrap around the GC. These 
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streams, known as the leading and trailing tidal tails, are made of matter tidally pulled 
away from the Sgr galaxy. They are on the orbital plane of the Sgr galaxy, which crosses the 
Galactic plane near the position of the Solar system. The leading tidal tail results from stars 
(and dark matter particles [see that were originally between the Sgr dwarf center and 
the center of the Milky Way. These stars that are closer to the GC orbit faster than the Sgr 
dwarf center, they move at lower gravitational potential and have shorter orbital periods 
around the GC. Thus, these faster orbiting stars (and dark matter particles) stream out 
ahead of the main body of the Sgr dwarf and create the leading tail. The stars that give rise 
to the leading tail are stripped at roughly a distance of 10 kpc from the GC, on the opposite 
side of the GC from the Sun. This distance of 10 kpc is determined by subtracting the 
current distance between the GC and the center of the Sgr main body, ~ 16 kpc, minus 
the Sgr tidal radius, ~ 6 kpc. This distance of 10 kpc determines the perigalacticon of the 
leading tidal tail, both where it was removed from the Sgr main body and again on this 
side (the solar side) of the GC. Similarly, the trailing tail is formed from stars that were on 
the opposite side of the dwarf when they were stripped, they are further out and thus move 
slower than the Sgr galaxy main body. 

The trailing stream roughly extends in a doughnut shape in the Southern galactic hemi- 
sphere around the position of the Sun, at a distance of about 20 kpc from us. The leading 
tail goes around the GC in the Northern galactic hemisphere, turns around at a radius of 
about 35 kpc, then approaches the galactic plane from the Northern galactic hemisphere 
and crosses this plane passing through the Solar system (see Fig. 11 of ^) in the general 
direction orthogonal to the Galactic plane. The streams carry dark matter as well as visible 
stars. Thus the leading stream brings additional dark matter into the vicinity of the Solar 
system. The effect of the leading Sgr stream for direct detection experiments, if the DM 
consists of weakly interacting massive particles, has recently been studied in ^2] where the 
estimated range of DM density in the stream was found to be 

Pstream = [2, 200] X lO^Mo/kpc^ = [0.001, 0.08] GeV/cm^ . (3) 

This density corresponds to a fraction x=(0.3-25)% of the local density of the isothermal 
Galactic halo for the usual estimate ph = 0.3 GeV/cm^, and to a fraction x=(0.45-37.5)% 
of the local halo density for a lower estimate ph = 0.2 GeV/cm^ (e.g. see the compilation 
of possible local halo density values for spherically symmetric halo models in [IE,). This 
evaluation assumes a constant density along the stream [T^ . 

The trailing stream may lead to some small increase of UHECR events coming from the 
Sgr dwarf plane, but this increase will be very small and would be more difficult to detect 
than the leading trail. Recall that the flux depends on the integral of the density over 
the line of sight (see Eq. (1)), and we can see along the axis of the leading trail but only 
transversally through the trailing tail. 
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We concentrate here on the leading tidal stream, which passes through the solar system. 
The width of the tidal stream in 2MASS M stars is estimated to be 4-8 kpc and the best-fit 
projection brings the center of the leading tail within 2 kpc of the Sun [H]. We can see in 
Fig. 11 of [0] that the stream starts turning around at about (25-35) kpc from the Sun, but 
for smaller distances can be approximated by a cylinder. For simplicity here we model the 
leading stream as a cylinder of radius -R=(2-4) kpc and height L= (25-35) kpc "above" the 
plane of the galaxy, i.e. on the side of the North Galactic Pole (NGP). We assume that the 
central axis of the cylinder passes through the position of the Sun. The NGP is located at 
5 = +27° 24' of declination and a =12^ 49™ of right ascension. 

We then compare the expected UHECR flux 0|| in a solid angle of aperture 9 in the 
direction of the NGP (i.e. along the axis of the cylinder), with the flux 0^ coming from an 
identical cone taken in a direction perpendicular to both the directions to the NGP and the 
GC (a cone with axis in the galactic plane). There are two directions with this property, 
one pointing to 6 = —48°, a =9^ and the other pointing to 6 = +48°, a =21^. We define 
the asymmetry amplitude 

AstreUG) = T^- (4) 

In a spherically symmetric halo without the Sgr stream there should be no difference in 
the fluxes from both directions, thus ^stream =1- The presence of the stream leads to an 
asymmetry ^stream >1- 

Taking the evaluation of the density of the leading stream from jl2] to be Pstream = XPh, 
where ph is the halo density at the location of the Sun, and assuming that the density is 
constant along the stream, we have computed v4stream(6'), using the halo model in Eq. (2). 
Notice that v4stream(6') depends on the local overdensity x due to the stream and not on the 
local halo density itself. In Fig. 1 we plot the column density as seen from the location of 
the Sun for the halo model in Eq. (2) with Rc= 5 kpc and for several of our leading stream 
models. The increase in column density in the directions close to the axis of the stream is 
clearly visible, except for the lowest values of x- 

In Figs. 2 and 3 we plot the asymmetry Astream{9) for different values of x and two values 
of i?c = 5 kpc and 10 kpc (Fig. 2 and 3, respectively). With R and L as the dimensions 
chosen for the cylinder representing the stream, Astrea.m{9) is constant for 9 < arctan(_R/L), 
and decreases rapidly with increasing values of 9, as can be seen in Figs. 2 and 3. The 
maximum asymmetry, Astrea.m{9 < 6.5°) = 1.7, is obtained for the largest cylinder, lowest 
halo density (i.e. highest x) and small core radius, namely for _R=4 kpc, L=35 kpc, x = 0.375 
and Rc=5 kpc. A larger core radius decreases the asymmetry to Astrea.m{9 < 6.5°) = 1.57. 
Reducing L to 25 kpc one obtains v4stream(^^ < 9°) = 1.5 for Rc = 5 kpc (and 1.38 for 
Rc = 10 kpc). Reducing the radius of the cylinder to R = 3 kpc and reducing x to 
X = 0.25, with L=25 kpc the asymmetry amplitude becomes Astrea.m{9 < 7°) = 1.36 (1.28) for 
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Rc = 5 (10) kpc. These are some of the largest expected asymmetries, while the asymmetries 
for the smallest values of x a^^e negligible. For example, ^stream < 1-005 for _R=3 kpc, 
L=25 kpc and x = 0-003. 

In order to detect the asymmetry Astream(6') above the three-sigma level we need a min- 
imum number of events within the two cones of aperture 6 that we are comparing, 

^^ = v^%S^>3. (5) 

^N\\ + V^tream + 1 

which means 

N > 9 (^^trcam + 1) 
(^stream 1) 

For the maximum value of the asymmetry, ^stream = 1-7, we need at least A^=50 events; 
for ^stream = 1-4, we need N= 135 events; for ^stream = 1-2, we need A^= 500 events; for 
^stream = 1-1, wc uccd at Icast N= fcw 10^ cveuts; and for ^stream < 1-01, we would need > 
a few 10^ events. These are the number of events within the cones considered. While the 
number of events increases as the solid angle of the cones increases, the stream asymmetry 
decreases for cones of aperture larger than arctan(i?/L). Thus there is an optimum cone 
aperture which maximizes the reach of each observatory. 

Let us see how large an overdensity x the Pierre Auger Observatory, EUSO and OWL 
would be able to measure. The reach of each observatory is shown in Fig. 2 and Fig. 3 with 
thicker solid (red) lines. As can be seen in the figures, the integration cones for which the 
stream is best detectable through the asymmetry have an aperture 6 of about 15° to 20°. 
For larger cones, although the minimum asymmetries that can be detected decrease with 
increasing number of events A^ within each cone, the predicted stream asymmetry decreases 
faster, so that the reach in x does not improve. 

To estimate the number of expected events, we need a value for the UHECR flux. At 
energies of 6 x 10^^ eV the AGASA and HiRes cosmic ray fluxes differ by about a factor of 
3 (see e.g. Fig. 1 in ^Hl)- Since it is the AGASA data that do not show a GZK cutoff, we 
adopt the AGASA value for the flux, with an energy dependence of which is typical of 
decaying dark matter models jB]. This gives us an isotropic UHECR flux integrated above 
6 X 10^^ eV of about 0.07 (km^ sr yr)~^. At higher energies the data points of AGASA, which 
are about one order of magnitude higher than those of HiRes, give an integrated isotropic 
flux of about 0.04 (km^ sr yr)-^ for E > 1 x 10^° eV. 

The southern station of the Pierre Auger Observatory (see for example ^M,) is already 
under construction since the year 2000 in the Province of Mendoza, Argentina, at a latitude of 
about 35.2° South. It will be completed by 2004. The northern station is proposed for Utah 
and would then enable continuous full-sky coverage. The southern array consists of 1600 
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particle detectors spread over 3000 km^ with fluorescence telescopes placed on the boundaries 
of the surface array. The expected angular resolution is less than 1° for all energies (0.5° 
at lO^^eV). This angular resolution is very good to localize the flux from the stream. The 
limiting effective aperture for the full southern array is 7350 km^sr, including only zenith 
angles less than 60°, which is the region of the sky where the reconstruction algorithm 
works best. Given that the latitude of the southern array is 35.2° South, this means that 
only objects at declination smaller than +24.8° can be reconstructed well. However, the 
stream is centered at the NGP, which has a declination of +27.4°, and so it is outside the 
best detection region of Auger South. 

Thus, the stream asymmetry could not be studied by the Southern Auger site unless 
events below a 60° zenith angle can be analyzed effectively. In this case the aperture would 
increase by about 50% ^3], becoming 11250 km^sr, and there would be a total of 790 events 
per year with E > 6x 10^^ eV. In the extreme case in which zenith angles as large as 85° could 
be reached, and we are not sure that this can be done, the number of events within the cones 
of solid angle n we are studying would be a fraction / = fi/(3.647r) = (1 — cos6')/(l + cos40°) 
of the total number of events observed, and only cones with aperture 6 < 22.4° would be 
completely included in the observable portion of the sky. In this case the minimum v4stream(^) 
observable by Auger South in 10 years as function of 6 for 6 < 22.4° can be seen in Figs. 2 
and 3 (thicker solid (red) line with the label Auger S). In particular, for 6 = 20°, we have 
/ = 0.034, and after 10 years we would get N = 268 events with E > 6 x 10^^ eV. This 
number of events would allow us to detect an asymmetry Astream(20°) > 1.28. From Figs. 4 
and 5, which plot the overdensity x due to the stream as a function of the asymmetry within 
a cone of 20° aperture, we see that if Auger South could observe at zenith angles of up to 
85°, it could detect the stream if x is at least x = 0.3 (this value occurs for i? = 4 kpc, L 
= 35 kpc and Rc = 5 kpc). 

For Figs. 4 and 5 we have chosen a cone aperture of 20° because this aperture is close to 
that one that maximizes the reach, although the latter depends somewhat on the observatory 
(see Figs. 2 and 3). 

If the Northern Auger site would become a reality, the NGP and thus the stream would 
be in full view of this new observatory. The projected total effective area of 44000 km^sr 
of the combined Southern and Northern Auger sites, with total Att sr sky coverage, would 
yield 3080 events with E > 6 x 10^^ eV. The number of events within one of the cones of 
solid angle fl we are studying is a fraction / = Q/{47t) = (1 — cos^)/2. For 6 = 20°, this 
gives / = 0.030. Thus in 10 years of combined Southern and Northern Auger data we would 
obtain = 924 events with E > 6x 10^^ eV. This would allow the detection of asymmetries 
above the thicker solid (red) line labeled Auger S+N shown in Figs. 2 and 3. In particular 
we could detect Astream(20°) > 1.14, which would translate into the detection of a stream 
with X > 0.14 (see Figs. 4 and 5). 
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The Extreme Universe Space Observatory (EUSO) jTH]; a mission on-board the Interna- 
tional Space Station, will image the air fluorescence of the extensive air showers produced 
by UHECR in the atmosphere. Its angular resolution is expected to be less than 2° and 
a conservative estimate of the duty cycle is 10%. It will operate for three years, with an 
effective trigger aperture which increases very fast with energy from 1.0 x 10^ km^ sr at 
E = A X 10^^ eV to 4.5 x 10^ km^ sr at E > 1 x 10^° eV (these effective apertures include 
the 10% duty cycle). The sky coverage is expected to be very close to complete, thanks to a 
suitable inclination of the EUSO detector axis with respect to the vertical and to the char- 
acteristics of the International Space Station orbit (the ISS orbit has a typical inclination of 
about 52° with respect to the ecliptic, a period of about a 16th of a day, and a precession 
of the line of nodes of about 5° per day ^7j). Thus, within An sr of sky coverage, we expect 
about 1800 events per year with E > 1 x 10^° eV. The number of events within the cones 
of sohd angle fl we are studying is a fraction / = Q/^An) = (1 — cos^)/2 (/ = 0.030 for 
6 = 20°). This means that in three years EUSO could detect asymmetries above the thicker 
solid (red) line labeled EUSO shown in Figs. 2 and 3. In particular, in three years of oper- 
ation = 162 events would be collected with E >lx 10^° eV for 6 = 20°. Thus, EUSO 
could measure a minimum asymmetry ^stream(20°) = 1.36, which is larger than the largest 
we obtained in our models at 6' = 20° (see Figs. 2 and 3). 

The Orbiting Wide-angle Light-collectors (OWL) ^H] mission would consist of a pair 
of satellites observing extended air showers in the atmosphere produced by UHECR with 
stereoscopic view. The angular resolution would be less than 1°. It would have an instan- 
taneous aperture rapidly increasing with energy, from 4 x 10^ km^ sr at -E = 3 x 10^^ eV 
to about 2 X 10^ km^ sr at E > 1 x 10^° eV. Taking into account that OWL can only view 
the dark side of the Earth, and accounting for the effect of the Moon, clouds and man- 
made light, the duty cycle is about 10%. In fact the continuous effective aperture would 
be 2.3 X 10^ km^ sr for > 1 x 10^° eV [TBj and about two orders of magnitude smaller 
at E = 3 X 10^^ eV. Because of the steep rise of the aperture with energy, OWL would 
see more events at E > 1 x 10^° eV than at E > 3 x 10^^ eV. With a flux of 0.04 (km^ 
sr s)~^, OWL would collect 9200 events per year. The number of events within the cones 
of solid angle fl we are studying is a fraction / = f2/(47r) = (1 — cos^)/2 of these 9200 
events per year. In five years of operation the minimum asymmetry OWL could detect is 
shown as a function of 6 in Figs. 2 and 3 with the thicker solid (red) line labeled OWL. 
In particular, OWL would detect 276 events per year in a cone of aperture 6 = 20°. This 
corresponds to a detection of = 1380 events within the same cone in 5 years. Thus OWL 
could detect stream asymmetries above 1.12, i.e. could reach local overdensities due to the 
stream x > 0.12. In terms of overdensity due to the stream, OWL's five-year reach is better 
than the reach of the combined Auger South and North in 10 years. 

In conclusion, EUSO would not be able to see the stream. The Southern Auger site alone 
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cannot study the stream asymmetry unless events below a 60° zenith angle can be analyzed 
effectively. In the hypothetical case in which zenith angles as large as 85° could be reached, 
in 10 years Auger South alone could test local overdensities in the stream of x > 0.3 (if 
Rc = 5 kpc). This is possible only if ph < 0.3 GeV/cm^. For example, if ph = 0.3 GeV/cm^, 
the minimum stream density testable by Auger South alone in 10 years would be Pstream = 
0.09 GeV/cm^, while if ph = 0.2 GeV/cm^, it would be p,tream = 0.06 GeV/cm^. 

The best chance of detecting the Sgr leading tidal stream in UHECR would come ei- 
ther with 5 years of observation with OWL or 10 years of observation with the combined 
Northern and Southern Auger sites. In fact, the combined Northern and Southern Auger 
data, integrated for about 10 years, would allow us to detect local stream overdensities of 
X > 0.14 (if Rc = 5 kpc). This means pstream > 0.042 GeV/cm^ if p^ = 0.3 GeV/cm^, or 
Pstream > 0.028 GcV/cm^ if p^ = 0.2 GcV/cm^. 

The OWL data, integrated for about 5 years, would allow us to detect local stream 
overdensitites of x > 0.12 (if Rc = 5 kpc). This means Pstream > 0.036 GeV/cm^ if ph = 
0.3 GeV/cm^, or pstream > 0.024 GeV/cm^ if ph = 0.2 GeV/cm^. 

The numbers given here are only indicative and we do not attempt to provide an estimate 
of their error. Our conclusions depend on the assumed UHECR fluxes which we took here 
to be at the level measured by AGASA. Clearly, the searches would become less sensitive 
for smaller fluxes. Our conclusions also depend on the halo model, in that the asymmetry 
due to the stream would decrease if the halo is flattened in the direction of the stream. 

P. G. would like to thank Benjamin Stokes for helpful information on the sky coverage of 
UHECR experiments. We thank the referee for suggesting corrections and Dmitry Semikoz 
for usefull discussions. The work of G .G. was partially supported by DOE grant DE-FG03- 
91ER40662 and by NASA grant NAG5-13399. 
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Figure 1: Halo and stream column densities for several stream models, for halo core radius 
Rc = 5 kpc, as a function of the angle ip measured from the Galactic Center in the direction 
of the North Galactic Pole (the Galactic Center is at = 0, the North Galactic Pole is at 
■0 = 90°, and the Galactic Anticenter is at = 180°). 
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Figure 2: Asymmetry amplitudes Astream(^) as a function of the cone aperture 9, for 
Rc=5 kpc and several values of the overdensity x due to the Sgr stream and the radius 
R and length L of the cylinder that models the Sgr stream. The thicker solid (red) super- 
imposed lines indicate the reach of each of the observatories. 
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Figure 3: Same as Fig. 2 but for Rc=lO kpc. 
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Figure 4: Local overdcnsity due to the stream x as a function of the asymmetry amphtude 
for a cone of 20° aperture, ^stream(20°), for i?c=5 kpc. 
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Figure 5: Same as Fig. 4 but for Rc— 10 kpc. 
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